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ABSTRACT 

We present a high spatial resolution (« 20 pc) of 12 CO(2 — 1) observations of the lenticular galaxy 
NGC4526. We identify 103 resolved Giant Molecular Clouds (GMCs) and measure their properties: 
size R, velocity dispersion a v , and luminosity L. This is the first GMC catalog of an early-type galaxy. 
We find that the GMC population in NGC4526 is gravitationally bound, with a virial parameter a ~ 1. 
The mass distribution, dN/dM oc M ~ 2 - 39±003 , is steeper than that for GMCs in the inner Milky Way, 
but comparable to that found in some late-type galaxies. We find no size-linewidth correlation for the 
NGC4526 clouds, in contradiction to the expectation from Larson’s relation. In general, the GMCs 
in NGC4526 are more luminous, denser, and have a higher velocity dispersion than equal size GMCs 
in the Milky Way and other galaxies in the Local Group. These may be due to higher interstellar 
radiation field than in the Milky Way disk and weaker external pressure than in the Galactic center. 
In addition, a kinematic measurement of cloud rotation shows that the rotation is driven by the 
galactic shear. For the vast majority of the clouds, the rotational energy is less than the turbulent 
and gravitational energy, while the four innermost clouds are unbound and will likely be torn apart by 
the strong shear at the galactic center. We combine our data with the archival data of other galaxies 
to show that the surface density £ of GMCs is not approximately constant as previously believed, 
but varies by ~ 3 orders of magnitude. We also show that the size and velocity dispersion of GMC 
population across galaxies are related to the surface density, as expected from the gravitational and 
pressure equilibrium, i.e. a v R~ x / 2 oc E 1 / 2 . 

Subject headings: ISM: clouds - galaxies: individual: NGC4526 - galaxies: ISM - galaxies: lenticular 
- radio lines: ISM. 


1. INTRODUCTION 

Giant molecular clouds (GMCs) are the sites of star 
formation in galaxies. The existing correlation be- 
twee n molecular gas surface density and star formation 
rate llWong fe Blitzll2002t iBigiel et al.1 120081 iLerov et~ahl 
I2013T) implies that the formation and evolution of GMCs 
are essential to understand the build-up of stellar masses 
in galaxies. However, up-to-date studies of extra- 
galactic GMC p opulations are lim it ed to Local Grou p 


galaxies; LMC (llFukui et al.l 2008t 

Wong et al. 

2011 ), 

SMC (Mizuno et al. 20017 M31 ( 

Rosolowsky 

2007), 

M33 (lEngargiola et al. 20031; Rosolowsky et al. 

20071 

Gratier et al.l 20121). and IC10 (ILerov et al. 20061). the 


nearby spirals; M64 dRosolowskv fc Blitdl2005H and M51 
( Colombo et al.l 1201 4l ) . and the nearby starburst; M82 
( Keto et al.ll2005fl and NGC253 (ILerov et alJl2014l) . due 


to the limited angular resolution and sensitivity of ra¬ 
dio telescopes. Galaxies in the Local Group are mostly 
dwarfs with few spirals. Therefore, additional study 
of GMCs in the early-type galaxies, such as NGC4526, 
is needed to provide a comprehensive analysis of GMC 
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properties across different galaxy environments. 

There are three resolved GMC properties that we can 
directly measure: size, linewidtli, and luminosity. The 
relationships b etwee n these properties were first stud¬ 
ied bv iLarsonl •: 1 98 H i. who suggested the importance of 
turbulence in the stability of GMCs agai nst s elf-gravity. 
These relations were then refined by Solomon et al. 
(119871 hereafter S87) for GMCs in the Milky Way disk. 
Basically, GMC properties in the Milky Way can be de¬ 
scribed by three Larson’s ‘laws 1 : (1) GMCs are gravita¬ 
tionally bound objects, (2) the size and velocity disper¬ 
sion of GMCs follows a a v oc R 0,5 relation, and by impli¬ 
cation (3) the surface density of GMCs is approximately 
constant (Egmc ~ 170 M 0 pc -2 , S87). Interestingly, 
these relations also hold true, albeit with sca tter, f or ex- 
tragalactic GMCs in Local Group galaxies dBlitz et al.1 
l2007t iBolatto et al.ll2008t iFukui fc Kawamurall2010H . At 

face value, this suggests that GMC properties are uni¬ 
versal. 

However, further studies reveal that G MC proper¬ 
ties c an deviate from Larson’s relations. iHever et all 
(120011 1 found that low mass GMCs (M < 1O 3 M 0 ) in 
the outer part of the Milky Way are not gravitation¬ 
ally bound. Their luminous masses, inferred from the 
CO-to-H 2 conversion factor, are smaller than their virial 
masses. The required external pressure to bind these 
clouds is Pext/k ~ 1 x 10 4 K cm” 3 . Furthermore, 
when re-examinin g S8 7 clouds using more sensitive in¬ 
struments, IHever et~aI1 ( 2009 ) found that the surface 
density of Milky Way disk GMCs varies from ~ 10 to 
200 Mq pc -2 , and they deviate from gravitational equi- 
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librium. iField et all (l2011i) suggest that these clouds 
may be in pressure virial equilibrium, where the clouds’ 
mass and radius are described by a Bonnor- Ebert sp here 
with various external pressures ( Bonnorl Il956t febertl 
195511. In addition, several authors (e.g. iKesrell Il989t 


Ballesteros-Paredes fc Mac Lowll2002fl argue that the ob¬ 


served constancy of surface density might be affected by 
observational biases. 

If this is really the case, then what factors set 
the different properties of GMCs? These parameters 
may be external (environ me ntal) effects, su ch as hy¬ 
drost a tic pressure (IElmegree.nl 1 199. 'll : iBlitz fc Rosolowskvl 
12004 iMeidt et al.1 I2013TK interstellar radiation field 
strength dMcKee 19891). and shear from g alaxy rotation 
ijKoda et all 120091 : iMivamoto et al.1120141 . or internal, 
such as feedback of the star formation that is embedded 
inside GMCs ( McKee! 119891 ). To answer this question, 
we need a complete sample of GMCs across different en¬ 
vironments: from bulge to spiral-arm and inter-arm re¬ 
gions, from late-type to early-type galaxies, and from low 
to high metallicity galaxies. 

In this respect, we analyze the GMC properties in 
the bulge of NGC4526, an SO-type galaxy in the Virgo 
cluster. NGC4526 is unusual because all of the CO in 
the galaxy has been observed at a linear resolution of 
ss 20 pc, sufficient to resolve Milky Way sized GMCs. 
The galaxy has pro mine nt cen tral dust lanes with mass 
~ 1O'M 0 (|Ciesla et all 120141) and supers olar metallic¬ 
ity flog {Z/Zq) ~ 0.2: IDavis et al.l l2013al) . but lack of 
star formation (SFR ss 0.03 M 0 yr” 1 ; lAmblard et al.1 
2014 and devoid of a tomic gas (Mhi < 1.9 x 1O 7 M 0 


Lucero fe Y oung 2013J). The HI deficiency in this galaxy 


may be caused by ram pressure or evaporation by hot gas 
(as the galaxy resides in the Virgo cluster), or by abrupt 
conversion of HI into molecular gas due to high pres¬ 
sure fElmegreen 1199,4 . There is no indication of recent 


tidal interaction (lYoung et al.ll2008l) , suggesting that this 
mechanism is not the primary cause of HI deficiency in 
NGC4526. 


In fact, the molecular gas in NGC4526 is confined 
within th e cent ra l ~ 1 kpc region (the top panel of Figure 
1: IDavis et al.l!2?H3b ). The central regions of galaxies, 
such as in the Milky Way, tend to have high interstel¬ 
lar pressures (~5x 10 6 K cm” 3 ; Snergel &; Blitzlll992D 
strong magnetic fields (~ 1 mG; lYusef-Zadeh fc Morrid 
179871 ). and lower than expected star formation rate 
(iLongmore et al.ll2013t ). These properties offer a unique 
environment for GMCs in NGC4526, significantly differ¬ 
ent than those studied in other galaxies. 

This paper is organized as follows. In §2 we de¬ 
scribe the data and methodology to identify GMCs in 
NGC4526. The properties and kinematics of the GMCs 
are reported in §3 and §4, respectively, and catalogued 
in Table 1. We discuss the pressure balance of GMCs in 
§5.1, the Larson’s ‘laws 1 in §5.2 and §5.3, and the envi¬ 
ronmental effects to the GMC properties in §5.4. Lastly, 
we summarize our findings in §6. 


2. DATA AND METHODOLOGY 
2.1. Data Descriptions 

NGC4526 was observed in the 12 CO(2 — 1) line (230 
GHz or 1.3 mm) using the Combined Array for Research 
in Millimeter-wave Astronomy (CARMA) in A, B, and C 


configurations (jBock et al.1 120061) . The data were taken 
as part of the mm-Wave Interferometric Survey of Dark 
Object Masses (WISDOM) pr oject. Resu l ts for t he in¬ 
nermost CO were reported by IDavis et~aD (j2013b ), who 
showed that the kinematics of the central CO imply the 
presence of a 4.5 x 1O 8 M 0 supermassive black hole. 

The beam width of the observations is 0.278 x 
0.173 arcsec 2 and the spectral resolution (after Han¬ 
ning smoothing) is 10 km s . This beam width covers 
5.56 x 3.46 pixels and corresponds to a projected physi¬ 
cal si ze of sa 22 x 14 pc 2 at the adopted distance of 16.4 
Mpc dTonrv et al.ll2001l) . These high resolution data en¬ 
able us to resolve individual GMCs and measure their 
properties, si nce the typ ical Milky Way’s GMC sizes are 
~ 50 pc (e.g. lBlit3il993t ). 

The noise in our data is not uniform, with higher noise 
appearing at the corners of the data cube. The overall 
pixel-by-pixel root-mean-square (rms) noise cr rms distri¬ 
bution is a positively-skewed Gaussian with minimum, 
average, and maximum values of 0.33, 0.71, and 1.33 K, 
respectively. 

2.2. Methodology 

We identify GMC candidates in NGC 4526 using the 
modified CLUMPFIND algorithm (Wil liams et al.1 
IT99i). implemented in the CPROPS program 
( Rosolowskv fe Lerovl I2006L hereafter RL06). The 
main goal of this program is to identify all real clouds 
and minimize false detections due to noise fluctuations. 
Descriptions of the CPROPS program, together with our 
chosen values of the input parameters of the program, 
are given in Appendix A and B. 

As a result of the CPROPS analysis, 241 GMCs are 
identified in NGC4526, of which 103 are resolved. We 
assume all GMCs are real since the probability of false 
detections is very small (Appendix C). In the bottom 
panel of Figure 1, we show the integrated CO emis¬ 
sion of connected regions that have brightness temper¬ 
atures Tb > 2cr rms and have at least one pixel with 
Tb > 3tr rms . The locations of resolved and unresolved 
clouds are marked as blue and red ellipses, respectively. 
The peak S/N distribution of identified GMCs is shown 
as an inset. The mean peak S/N of resolved and unre¬ 
solved clouds is 5.6 and 4.7, respectively. 

Most GMCs are located in the molecular ring, a few 
hundred parsecs from the galactic center (Figure 1). This 
molecular ring is the largest contiguous CO emission in 
our data. In addition, there are a few clouds located 
in the central region. The outer region of the molecu¬ 
lar gas exhibits a spiral-arm structure, possibly with an 
outer pseudo-ring, which fragments into smaller struc¬ 
tures consisting of one or multiple GMCs. Except for 
one cloud that is described below, all identified GMCs 
are within 900 pc of the galactic center, i.e. inside the 
bulge of NGC4526. Note that our primary beam covers 
all of the CO emission in the galaxy, so our GMC catalog 
is complete. 

There is one unresolved cloud which is located on the 
edge of the data cube (not shown in Figure 1), out of the 
plane of the CO emission. We overplotted the location 
of this cloud with the HST archival image of the galaxy. 
Although this cloud is likely to be real (Appendix C), its 
distance is uncertain, i.e. it may be located outside the 
galaxy. We exclude this cloud from the following analysis 
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Fig. 1.— Top: The integrated intensity map of NGC4526. The map is created by applying a Gaussian fit to each spectrum in the data 
cube. We exclude any Gaussian that has peak less than 2.5<7 rms . The ellipses divide the CO emission into three zones: inner region, 
molecular ring, and outer region. Bottom: Identified GMCs in NGC4526 overploted on the masked integrated intensity map. The mask 
covers regions with connected emission above 2cr rms and having at least one pixel with 3cr rms . The beam size and the projected physical 
size are indicated. The blue and red circles mark the location of the resolved and unresolved clouds, respectively. The distribution of the 
peak S/N of clouds is shown as an inset in the bottom right corner. 


but keep it in the catalog (as cloud no. 80 in Table 1). 
Inclusion of this cloud does not alter the conclusions of 
our analysis. 

3. CLOUD PROPERTIES 
3.1. Definition of GMC Properties 

Cloud properties, such as position, size, velocity dis¬ 
persion, luminosity, and mass, are catalogued in Table 
1. Here, we briefly describe the method used to measure 


the cloud properties. Full explanations of the method 
are given in RL06. 

The cloud size, R , is measured using the deconvolved 
second-moment: 

R = V \! (<4aj[° K] - ^eam) 1/2 «m[° K ] ~ ^beam) 172 - 

(!) 

where 77 is a factor that depends on the density distribu¬ 
tion of spherically symmetric clouds. A uniform sphere 
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TABLE 1 

NGC4526 Cloud Properties 


ID 

RA(2000) 

[h:m:s] 

Dec(2000) 

[° 

Vlsr" 

R b 

[pc] 

SR b 

[pc] 

( T v a 

(5(T„ a 

L c 

SL C 

Mi um a 

<5Af lum a 

S/N 

-'Hi, max 

[K] 

^shear 

[pc] 

1 

12:34:3.5 

7:41:54.7 

276.0 

24.00 

7.65 

6.90 

0.48 

3.76 

1.24 

1.88 

0.62 

5.0 

7.7 

0.52 

667 

2 

12:34:3.5 

7:41:54.9 

283.4 



6.14 

1.30 

0.67 

0.60 

0.33 

0.30 

4.1 

6.3 

0.55 

619 

3 

12:34:3.4 

7:41:55.1 

286.7 



6.37 

1.21 

0.27 

0.19 

0.14 

0.09 

3.4 

4.7 

0.57 

587 

4 

12:34:3.6 

7:41:54.3 

279.6 

22.16 

5.54 

5.42 

0.55 

2.35 

0.70 

1.18 

0.35 

6.6 

10.4 

0.44 

798 

5 

12:34:3.5 

7:41:55.2 

286.7 



5.64 

1.55 

0.71 

0.68 

0.35 

0.34 

4.6 

6.7 

0.57 

583 

6 

12:34:3.4 

7:41:55.4 

293.2 



9.74 

2.35 

1.13 

0.60 

0.57 

0.30 

5.2 

6.7 

0.58 

563 

7 

12:34:3.5 

7:41:55.4 

298.3 

22.01 

5.52 

7.23 

0.96 

2.39 

0.76 

1.20 

0.38 

5.7 

8.1 

0.56 

599 

8 

12:34:3.4 

7:41:55.1 

298.5 



6.28 

1.67 

0.67 

0.99 

0.33 

0.50 

3.6 

5.1 

0.58 

559 

9 

12:34:3.4 

7:41:55.6 

297.2 

29.60 

6.81 

7.25 

0.63 

5.75 

1.35 

2.88 

0.68 

6.2 

9.0 

0.61 

492 

10 

12:34:3.3 

7:41:55.5 

311.0 

12.11 

9.12 

8.01 

1.41 

1.42 

1.35 

0.71 

0.68 

4.4 

6.4 

0.61 

480 

11 

12:34:3.3 

7:41:55.7 

314.9 



6.27 

0.55 

1.05 

0.68 

0.53 

0.34 

5.1 

7.1 

0.62 

440 

12 

12:34:3.3 

7:41:56.1 

309.0 

28.58 

5.57 

7.12 

0.33 

5.52 

1.16 

2.77 

0.58 

7.3 

10.2 

0.63 

424 

13 

12:34:3.3 

7:41:56.4 

323.3 

31.38 

5.24 

7.81 

0.86 

5.84 

1.08 

2.93 

0.54 

7.3 

10.3 

0.66 

348 

14 

12:34:3.3 

7:41:55.2 

323.5 

13.38 

8.55 

8.21 

4.00 

2.01 

0.87 

1.01 

0.44 

5.1 

7.4 

0.59 

543 

15 

12:34:3.3 

7:41:56.4 

333.7 

18.72 

5.54 

9.56 

2.32 

3.18 

0.87 

1.59 

0.44 

6.6 

8.9 

0.64 

388 

16 

12:34:3.2 

7:41:56.5 

344.0 

16.65 

7.01 

8.75 

2.44 

2.00 

0.83 

1.00 

0.41 

5.8 

7.8 

0.76 

293 

17 

12:34:3.2 

7:41:56.8 

328.9 



7.64 

1.41 

0.29 

0.21 

0.14 

0.11 

4.1 

5.8 

0.70 

320 

18 

12:34:3.4 

7:41:56.1 

335.7 



7.39 

3.21 

0.80 

0.91 

0.40 

0.46 

4.4 

6.0 

0.60 

499 

19 

12:34:3.3 

7:41:56.2 

330.3 

10.48 

7.41 

8.53 

1.50 

1.47 

0.64 

0.74 

0.32 

5.5 

7.9 

0.66 

360 

20 

12:34:3.2 

7:41:56.3 

337.5 

7.57 

11.12 

6.44 

1.74 

0.96 

0.97 

0.48 

0.49 

3.9 

5.9 

0.68 

336 


Note. — Table 1 is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance 
regarding its form and content. 
a Units are km s~ 1 . 

b The size of unresolved clouds is less than the linear size of the beam, denoted as fl = .... 
c Units are 10 5 K km s -1 pc 2 . 
d Units are 10®Mg. 
e Units are km s —1 pc — b 

f Distance from the center of NGC4526, assuming clouds are in the plane of the galaxy with axis ratio of 0.216 and position angle of 20.2°. 


has rj = y/5, while an isothermal sphere has rj = 3. Here, 
we adopt 77 = 1.19, a value from S87, to make it con¬ 
sistent with previous studies. The major, cr ma j, and mi¬ 
nor, crmin, dispersions are the spatial second moments, 
weighted by the intensity, along the major and minor 
axis of the clouds, respectively. Both <r ma j and cr m i n are 
extrapolated to zero intensity (0 K) to avoid bias due to 
the limited sensitivity of the instrument. Deconvolution 
is applied, by the inclusion of the dbeam terms in equa¬ 
tion ( 1 ), to avoid bias due to the finite beam resolution 
(RL06). Since the beam is not circular, we take hbeam 
as the geometrical mean of the major and minor axes 
of the beam. The uncertainty, SR, is determined using 
bootstrap resampling. 

In order to measure the clouds’ velocity dispersions, 
we attempt Gaussian fitting to the composite spectrum 
of each cloud through the following steps. First, we cal¬ 
culate the offset of the mean velocity at all lines of sight 
within the cloud ( Xi,yt ), with respect to the mean ve¬ 
locity at the center of the cloud (xo,j/o)- This offset is 
caused by large scale motions, such as the cloud’s rota¬ 
tion or shear due to galactic rotation. Then, we shift 
each line of sight velocity spectrum to match the mean 
velocity of the central position of the cloud. Except for 
the innermost clouds, this shift removes any velocity dis¬ 
persions due to large scale motions, and leaves only tur¬ 
bulence and thermal broadening as sources of velocity 
dispersions. To make a composite spectrum, we take 
the average velocity profile from each line of sight. Fi¬ 
nally, we fit the composite spectrum with a Gaussian. 
The standard deviation of the Gaussian fit is taken as 


the velocity dispersion, a v , of the cloud. We take the 
uncertainty, Sa v , from the bootstrap resampling. 

R.L06 shows that measurements of velocity dispersion 
suffer bias towards higher values due to the finite spectral 
resolution of the instruments. Therefore, we take the 
deconvolved value, cr t ,,d c , to remove this bias using the 
same prescription as in RL06: 


&v,dc 



( 2 ) 


where Sv is the spectral resolution, and Sv (27t) —0 5 is 
the standard deviation of a Gaussian that has an inte¬ 
grated area equal to a spectral channel with width Sv. 
For brevity, we refer to cr Vt dc as a v . Note that the de¬ 
convolved values of a v and R are always smaller than 
their measured values (c.f. equation 1 and 2). Thus, any 
clouds that are barely resolved would have the decon¬ 
volved value smaller than the resolution. 

The cloud luminosity is the integrated CO flux over the 
position-position-velocity volume occupied by the cloud 
(RL06): 

Lew-.) = ( 18Q * 3600 ) D 2 !>,. fa Sy Sv. (3) 


where Tb,* is the brightness temperature of the ith pixel 
in K, Sx and Sy are the pixel sizes in arcsec, D is the 
distance to NGC4526 in pc, and Lq.o is the cloud lumi¬ 
nosity in K km s ^ 1 pc 2 . The luminosity is extrapolated 
to 0 K intensity as described in RL06. The uncertainty, 
5L, is determined using a bootstrap resampling method. 
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log (M lum /M 0 ) 


Fig. 2.— Correlation between virial and luminous masses of 
GMCs. The solid line is the one-to-one relationship. The distribu¬ 
tion of log a with a log-normal fit is shown as an inset. The mean 
of the log-normal fit is a « 0.99 with a standard deviation of 0.14 
dex. Thus, the GMC population in NGC4526 is in gravitational 
equilibrium. 


There is an additional ~ 20% (absolute, but systematic) 
flux calibration uncertainty that we do not include in the 
analysis. 

The luminosity is then converted to mass using the 
Milky Way’s CO(l — 0)-to-H 2 conversion factor Xco = 


2 x 10 20 cm'" (K km s _ ) , which is assumed to be 

constant throughout the galaxy. This is a reason¬ 
able assumption since Xco does not vary significantly 
due to met allicity in the supersolar metallicity regime 
dBolatto et al.l [20131 and references therein). The ratio 
of 12 CO fl - 0) to 12 CO(2 - 1) intensity in NGC4526 
is 1.15 ( Crocker et all (2012 5. We refer to this Xco- 
derived mass as the luminous mass: M\ um /M Q = (4.4 x 


1.15) (Tco( 2 -i)/K km s^ 1 pc 2 ), which takes into account 
the mass contribution by helium. The factor of 4.4 comes 
from our adopted Xco value. 

We calculate each cloud’s distance from the center of 
NGC4526 using the assumption that they are located 
in the plane of the galaxy with axis ratio of 0.216 and 
position angle of 290.2°. The position angle is measured 
from the north, counter-clockwise to the receding part of 
the kinematical major axis of the galaxy. The axis ratio 
and position angle are calculated using p art of the multi- 
Gaussian expansion (MGE) fit of iGappellaril (120021 ). 


3.2. Gravitational Equilibrium of Clouds 

The mass of a gravitat ionally bound cloud is given by 
dBertoldi fc McKeelll992 . hereafter BM92) 


5<r 2 i? 

G 


1048 Mq 




(4) 


where a v is the 1-D velocity dispersion of the CO line. 
We refer to this mass as the virial mass. If the lumi¬ 
nous mass is equal, or comparable, to the virial mass, 
then the cloud is in gravitational equilibrium, where the 
kinetic energy balances its self-gravity. If the luminous 
mass is smaller than the virial mass, then, in addition 
to gravity, the clouds must be held together by the ex¬ 
ternal pressure of the ambient medium, P ext , to reach 
dynamical equilibrium. Such clouds are pressure-bound 
clouds. 

We can also define the virial parameter as the ratio 
between twice the kinetic energy and the gravitational 
energy, 

5 a 2 v R = M W 
GM\ um Mi um 

According to BM92, clouds with a ~ 1.13 are 

gravitationally-bound and clouds with a 1 are 
pressure-bound. In Figure 2, we plot the luminous ver¬ 
sus virial masses of the resolved clouds, together with the 
distribution of log a as an inset. A log-normal fit to the 
distribution yields a mean a = 0.99 ± 0.02 and a stan¬ 
dard deviation of 0.14 dex. Roughly 99% of the resolved 
clouds have 0.3 < a < 3, and 89% of the resolved clouds 
have 0.5 < a < 2. Therefore, the GMC population in 
NGC4526 is in a state of gravitational equilibrium. 

In equation (4), we assume that all clouds are spher¬ 
ically symmetric and have a uniform density distribu¬ 
tion. If the clouds were isothermal spheres, then the 
virial masses would be 60% lower than our calculation. 
Moreover, the uncertainty in Xco also affects the lu¬ 
minous mass measurements. From observations of local 
galaxies, the typical uncertainty in Xco is about 0.3 dex 
dBolatto et al.lI20T3L and references therein). Taken all 
together, this introduces an uncertainty of ~ 0.5 dex in 
the worst case. 

Variations of the input parameters of the CPROPS 
program do not affect our result that the cloud popula¬ 
tion is gravitationally-bound. However, we have to keep 
in mind that we do not yet take into account the mag¬ 
netic pressure and rotation of the clouds. The effect of 
rotation is discussed in §4.3. Furthermore, we find no 
correlation between the mass of a GMC and its distance 
from the galactic center, possibly because all GMCs are 
distributed in a small region (within a radius of 900 pc) 
inside the bulge of the galaxy, so environmental varia¬ 
tions from GMC to GMC, such as ambient pressure and 
the interstellar radiation field, are small. 

3.3. Cloud Mass Distribution 

We use the luminous mass to determine the mass func¬ 
tion because it is well-defined even for unresolved clouds. 
Since the GMC population in NGC4526 is in gravita¬ 
tional equilibrium (Mi um ss M v j r ), we should not expect 
variation of the mass function between the two mass mea¬ 
surements. We determine the mass function using three 
different methods: equal bin-width, equal bin-size, and 
the cumulative distribution function. All measurements 
are taken from the most massive clouds in our sample 
(Mi um ~ 5.9 x 10 6 Mq) down to the completeness level 
of the observations. 

To determine the completeness level of our observa¬ 
tions, we create simulated Gaussian clouds. Their prop¬ 
erties are related through known scaling relationships: 
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log (M/M 0 ) 

Fig. 3. — Fraction of recovered clouds as a function of cloud 
mass, from the simulation described in the text. Simulated clouds 
with M > 5 x 10 5 Mq are well recovered by the program (solid 
line), contaminated by only a small fraction of false clouds (dashed 
line). The fraction of false clouds is negligible for low mass regime 
because these clouds are too small and faint, and hence undetected 
by the CPROPS program. We adopt log(M/MQ) = 5.7 as the 
completeness level of our observations. 

a v ex R 0 5 and M ex R 2 (e.g. S87). To mimic the ob¬ 
served data, we add the typical noise of our observations 
into the simulated data cubes. In total we consider 1600 
mock clouds with log (M/Mq) ranging from 4.9 to 6.7, 
with an increment of 0.2, and feed these mock clouds into 
the CPROPS program. A cloud is defined as recovered 
if its location in the data cube is within one beamwidth 
of its input location. Otherwise, this cloud is defined as 
false detection. The false detection rate is effectively zero 
for the least massive clouds because these false clouds, if 
they are exist, are too small to be recovered. We find 
that clouds with log(M/M 0 ) > 5.9 are well recovered 
by the program, while more than 80% of clouds with 
log (M/Mq) ps 5.7 are recovered (Figure 3). Therefore, 
we adopt the completeness level as log (M/Mq) =5.7. 


TABLE 2 

Best-fit parameters of the cumulative mass distributions 


Region 

Distance [pc] 

X 

M 0 [10° M 0 ] 

N 0 

All 

0 < d < 900 

-2.39 ±0.03 

4.12 ± 0.08 

9.40 ± 0.70 

Inner 

0 < d < 170 

-1.40 ±1.19 

1.88 ± 0.23 

9.13 ±33.21 

Ring 

170 < d < 580 

-2.38 ±0.03 

4.66 ±0.11 

5.66 ± 0.47 

Outer 

580 < d < 900 

-2.46 ±0.12 

2.56 ±0.12 

4.67 ± 1.12 


From these variations of x, we conclude that the equal 
bin-width method has large scatter due to the choice of 
the number of bins. 


3.3.2. Equal bin-size 

iD’Agostino fe Stephensl (|1986i ) and 

iMafz Apellaniz fc Ubedal (12005H found that vari¬ 
able bin-widths with equally divided numbers of data 
points per bin can minimize the binning uncertainty, 
and hence, more robust. This is because no bin has a 
much smaller number of data points than the others, 
in contrast to the equal bin-width method. For our 
data, the actual number of data points in a bin is not 
exactly the same; it can differ from that in other bins 
by one data point, due to non integer numbers after 
division. To check the robustness of the result, we vary 
the number of bins as in the equal bin-width method, 
and fit the resulting histogram with a straight line. 
The results are indeed more robust than the equal 
bin-width method, with a maximum slope of —2.31 and 
a minimum slope of —2.41 (left panel of Figure 4). 


3.3.3. Cumulative distribution 

In addition, we also calculate the mass function us¬ 
ing the (truncated) cumu l ative distribution function (e.g. 
iWilliams fe McKeg|1997t lRosolowskvll2005lf . 


N(M' > M)=N 0 



( 6 ) 


3.3.1. Equal bin-width 

In the equal bin-width method, we group the masses 
into bins of equal width in log-space. Then, each his¬ 
togram is fit with a straight line, weighted by the uncer¬ 
tainty of dN/dM , from the highest mass bin down to the 
completeness level. The slope of the best fit line, x, is 
the exponent of the mass function dN/dM oc M x . 

We calculate the uncertainty on the number of clouds 
in each bin as follows. First, the uncertainties of the 
masses SM are calculated through a bootstrap resam¬ 
pling method (RL06). Then, we use SM to calculate 
the uncertainty in dN/dM using Monte Carlo simula¬ 
tions. In these simulations, we resample the masses of 
the cloud, given a log-normal probability function with a 
mean M and a standard deviation SM. The resampled 
masses are grouped into the same mass bins as the data, 
so that each simulation gives a new mass distribution. 
We repeat these steps 10,000 times, and take the un¬ 
certainty in dN/dM as the standard deviation of these 
10,000 simulations. 

To check the robustness of our results, we vary the 
number of bins from 6 to 20. We find that x varies from 
—2.67 to —2.16 (left panel of Figure 4). The uncertainty 
of the slope is taken from the covariance matrix of the fit. 


and the simple (non-truncated) power-law distribution 
function, 

, f M\ X+1 

N ^ M > M) = (mo) ’ < 7 > 

where Mo is the cut-off mass of the distribution and Nq 
is the number of clouds with M > 2 1 ^ x+1 ' ) Mq, i.e. the 
cut-off point of the distribution. The truncated mass 
distribution takes into account that the mass distribu¬ 
tion of a population lack clouds more massive than Mq. 
The cumulative distribution function is robust against 
the number of bins since we do not bin the data into a 
histogram. 

We fit the cumulative distributions of masses us- 
ing the orthogonal dist ance regression method in Scipy 
dBoggs &; Rogersl I1990H . The fit is made for all data 
above the completeness level. We find x = —2.39 ± 0.03 
and Mo = (4.12 ± 0.08) x 10 6 Mq, in agreement with 
the equal bin-size method. The data are inconsistent 
with simple (non-truncated) power-law mass distribu¬ 
tions (middle panel of Figure 4). We find that the 
GMC mass distribution in N GC45 26 is steeper than 
in the inner MW (x = —1.5; iRosolowskvl l2Q05h . but 
comparable to the GMC mass distribution in the outer 
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Fig. 4.— Left: Best fit slope of the mass function against the number of bins. The equal bin-width method (blue squares) has large scatter 
due to the choice of the number of bins, while the equal bin-size method (red circles) is in agreement with the cumulative distribution mass 
function (dashed line). Middle: Fits to the cumulative mass function with and without truncation. The data favor the truncated (solid 
blue curve) over the non-truncated (dashed blue curve) mass function. Right: Cumulative mass distribution of the inner region (blue), 
molecular ring (green), and outer region (red), with the overlaid truncated fits. 


MW (x = —2.1; iRosolowskvl 12 0051 and central M33 
(x = — 2.0: iRosolowskv et al.ll2007blf . 

All three methods of measurements suggest that x < 
—2. In this case, most of the mass resides in the low mass 
clouds. Furthermore, the total mass diverges for integra¬ 
tion down to an infinitely small mass. Hence, there must 
be a lower limit to the cloud masses or a change in the 
slope, i.e. x > — 2 for lower mass clouds below our com¬ 
pleteness level, so that the total mass remains finite. 

We further divide the galaxy into three distinct regions: 
inner (0 < d < 170 pc), molecular ring (170 < d < 580 
pc), and outer (580 < d < 900 pc) region (concentric el¬ 
lipses in Figure 1), and measure their mass distributions. 
The molecular ring and the outer region have a similar 
mass distributions with an exponent of x ~ —2.4, while 
the inner region is much flatter (x ~ —1.4), albeit with 
a large uncertainty due to the small number of clouds in 
the inner region. The best-fit parameters are compiled 
in Table 2. A radial depe ndence o f the mass function 
was also discovered in M33 (iGratier et al.ll2012lf and M51 
dColombo et al.ll2014H . 

We calculate the total mass of detected GMCs (includ¬ 
ing the non-resolved clouds) to be Mqmc = (2.0 ± 0.1) x 
1O 8 M 0 . The total H 2 mass is M(H 2 ) = (3.8 ± 1.1) x 
1O 8 M 0 (lYouner et al.ll200M so the fraction of molecu¬ 
lar mass that residing in GMCs is Mqmc / 'M (H 2 ) ~ 0.53. 
This value is formally a lower limit, since there are GMCs 
with masses below the completeness level of our observa¬ 
tions that are undetected. The rest of the molecular gas 
may be in the form of diffuse gas which is undetected by 
interferometric observations. 

3.4. Larson’s Relations 

iLarsonl (1 9811 ) found that the velocity dispersion of 
GMCs is correlated with their size through a power-law 
relation with exponent of ~ 0.3. This correlation is sim¬ 
ilar to that expected if turbulence governs the velocity 
dispersion within clouds as described by the Kolmogorov 

6 We recalculate the total H 2 mass using Xqq as in this paper. 


law. In subsequent work, S87 refined the exponent to be 
0.5 ± 0.05 for GMCs in the Milky Way’s inner disk. 

Larson’s relations in the Mil ky Way consist of two in¬ 
dependent equations (e.g. S87. lBolatto et al1l2008H : 

/ D \ 

<j v « 0.72 (—- J km s -1 (8) 

and either 

/ 77 > \ 2.5 

Leo ~ 25 f —J K km s _1 pc 2 (9) 

or 

Leo ~ 130 f — (Tu _ 1 ~) K km s _1 pc 2 . (10) 
V km s J 


Since the linewidth, AV, is just AH = cr v \J 8 In 2, we 
refer to equation (8) as the size-linewidth relation. For 
extragalactic clouds in the Local Group, Bolatto et al. 
(2008) found a v oc R o e , Leo oc R 2 ' 54 , and Leo oc ex 3 ' 35 , 
which is close to the Milky Way relations. 

Interestingly, we find no size-linewidth correlation for 
NGC4526 (left panel of Figure 5), which is parameter¬ 
ized by very weak Pearson and Spearman correlation co¬ 
efficients (r ps = —0.18 and r s n = —0.14). This r esult 
is in l ine with GMCs in M 33 (|Gratier et al.l I2012H and 
M51 ( Colombo et al .) 201 1:. where no clear trend was 
observed (r sp = 0.12 and 0.16 for M33 and M51, respec¬ 
tively). The NGC4526 data are located above the Milky 
Way’s relation, which means that for a given size, GMCs 
in the bulge of NGC4526 have a higher velocity disper¬ 
sion than those in Milky Way disk GMCs by a factor of 
~ 3. This could be due to an environmental e ffect, since 
iShettv et al.l (|2012t l and lColombo et all (|2014tl found evi¬ 
dence that GMCs in the central regions of the Milky Way 
and M51 have a higher velocity dispersion than those in 
the disks. We discuss this environmental effect in §5.4. 

We also plot cloud luminosity against velocity disper¬ 
sion and size in the middle and right panels of Figure 5. 





















































































Fig. 5.— GMCs properties (radius, velocity dispersion, and luminosity) are plotted relative to one another. The color codes are for inner 
GMCs (distance < 65 pc; red squares) and outer GMCs (distance > 65 pc; blue dots). The fits of all data points (including the inner 
clouds) are shown as solid lines. The typical uncertainty is shown as a cross sign in the corner of each panel. For clarity, we do not plot 
the uncertainties of individual data. However, we fit the data points by including the non-uniform error bar of the individual data point, 
not only the typical uncertainty. The dashed lines are Larson’s relations for the Milky Way disk (S87), and the dotted lines are Larson’s 
relations with different normalization factors: 3, 0.03, and 5 from left to right panel, respectively. It shows that GMCs in NGC4526 are 
more turbulent and more luminous than equal-size clouds in the Milky Way disk. There is no size-linewidth relation, in contradiction to 
the expectation from Larson’s relation (left panel). 


The results of the error-weighted fit are 

Leo = 2258.5™ s ° (—) K km s' 1 pc 2 (11) 

and 

/ >. 3 . 2 ± 0.5 

Leo = 381.3^238 4 ( r^T ) K km S_1 P° 2 ’ 

V km s J 

( 12 ) 

which is shallower than the exponents in the Milky Way 
relations (equations 9 and 10). The correlation co¬ 
efficients are moderate for the luminosity-size relation 
(r ps = 0.63 and r sp = 0.67), and weak for the luminosity- 
velocity dispersion relation (7’ ps = 0.33 and r sp = 0.36). 

The quoted results above take into account all resolved 
clouds. In Figure 5, the inner clouds with distance < 65 
pc from the galactic center tends to have higher velocity 
dispersion, which may be due to contamination by the 
galactic shear (see §4 for detailed discussion of cloud kine¬ 
matics). The conclusions for the size-linewidth and size- 
luminosity relation are not affected if we exclude these 
inner clouds. However, the slope of linewidth-luminosity 
relation become steeper by excluding those inner clouds 

(Leo oc <t))- 8±0 - 6 ). 

In the middle panel of Figure 5, the clouds in NGC4526 
lie below the Milky Way luminosity-velocity dispersion 
relation. Thus, for a given velocity dispersion, the clouds 
are less luminous than GMCs in the disk of the Milky 
Way. Because CO luminosity is a tracer of the amount 
of molecular gas, the clouds in NGC4526 are more tur¬ 
bulent per unit mass than those in the Milky Way. Also, 
from the right panel of Figure 5, clouds in NGC4526 
are ~ 5 times more luminous than equal size clouds in 
the Milky Way, which means they have a higher surface 
density. Th is re sult is in agreement with GMCs in M51 
ijColombo et al.ll2014f) , where GMCs in the central region 
are brighter than those in the inter-arms region. Again, 
this could be due to environmental effects which are dis¬ 


cussed in §5.4. 

Finally, we find no bias that affects the results due to 
the choice of input parameters of the CPROPS program. 
Even though the properties of individual clouds vary by 
changing the input parameters, the overall distributions 
are similar (see Appendix D for details). 

4. CLOUD KINEMATICS 
4.1. Velocity Gradient of the Cloud 

Previous studies (e.g. Ilmara et all 120111: 
Rosolowskv et all 1 20 (Kit IPhillipsI fl999t iKane fe Clemensl 
1997t iCoodman et al .11 19931 ) have found velocity gradi¬ 
ents across atomic and molecular clouds and clumps, 
which are interpreted as rotation of the GMCs. More¬ 
over, most of the GMC rotation exhibits solid body 
rotation. Here, we perform analyses on the resolved 
clouds to gain insight into the origin of their kinematics 
and the role of rotation in the dynamical stability of 
GMCs. 

In order to quantify any rotation signature, we do the 
following. First, the velocity field (first-moment map) 
of the cloud is smoothed with a Gaussian kernel, where 
the dispersion of the Gaussian kernel is half the telescope 
beamwidth. The aim of this smoothing is to ‘average 1 the 
velocity field at the cost of losing independence among 
the neighboring pixels (left panels of Figure 6). Then, 
we fit the first-moment map of i ndividual clouds with a 
plane (e.g.lGoodman et a.l .11 1 9931 iRosolowskv et al .1120031 

Ilmara et all 120 1 lib 

vi os = v 0 + a(x - x 0 )+ b(y ~ y 0 ), (13) 

where (xo, yo) is the cloud’s central pixel coordinate, and 
dv dv 

a = IS “ d b= ai (14) 

are the velocity gradients along the x and y axes, vq, a, 
and b are free parameters to be determined from the fit. 
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Fig. 6. — Examples of plane fitting to find the rotation signature 
of GMCs. The right panels are GMC first-moment maps with the 
rotation axis (black line) overplotted. On the left panels, the mean 
velocity of each pixel is plotted against its perpendicular distance 
from the rotation axis. The black line is the fit and blue dots are 
the mean of the velocity in bins of perpendicular distance from the 
rotation axis (i.e. v± vs d). The cloud in the top panels shows 
the signature of solid body rotation with j. ~ 0.18, while the 
cloud in the bottom panels shows bow-shock motions and is well 
fitted by a parabolic curve with Xu,par a ~ 0-68- 

The angle from the positive a:-axis to the receding part 
of kinematical major axis of the cloud is tan -1 (6/a), and 
hence, the angle to the cloud rotation axis (i.e. the an¬ 
gular momentum vector) is 0 = tan -1 (6/a) + 90° (i.e., 
right hand rule). 

The angular speed of the cloud is given by 

n msw =\Z(£) 2+ (S) ’ <i5) 

where if is the angle from the cloud rotation axis to the 


sky plane. Since we can not measure if directly, we drop 
the cos ( if ) term from equation (15). Thus, the true an¬ 
gular speed is underestimated by a factor of cos (ip), i.e. 

^projected = f^true COs(^). 

The next step is to check whether the clouds show 
solid-body or differential rotation. We plot the mean 
velocity of each pixel within a cloud against its perpen¬ 
dicular distance from the cloud rotation axis, i.e. v± vs. 
d (Figure 6), and then we fit the data with a straight 
line. Solid body rotators should show a clear linear be¬ 
havior on this plot, where the constant slope is the an¬ 
gular speed, i.e. v± = ft d. On the other hand, the slope 
of Keplerian rotators varies with distance from the cen¬ 
ter, i.e. the slope gets shallower outside (resembles an 
S-shapc; v± oc cG 0,5 ), while bow-shock motion s tend to 
have a parabolic shape dKane Sz Clemenslll997f) . 

Finally, we divide the clouds into two groups: clouds 
that show solid-body rotation (SB) and clouds that show 
deviations from solid-body rotation (NSB). This division 
is based on the reduced y 2 value from the binned cj_ 
(blue squares in Figure 6) to the straight line fit. We 
classify clouds with y 2 lirle < 1.5 as SB and the rest as 
NSB. For NSB, we also fit v±(d) with a parabolic curve 
and calculate its y 2 para . Eye inspections confirm that 
this classification is reasonable. As a result, 46 of 103 re¬ 
solved clouds are classified as SB, while the rest are NSB. 
Examples of SB and NSB clouds are shown in Figure 6. 

4.2. Origin of Velocity Gradients 

A purely rotating galaxy with an inclination angle i 
has a line of sight velocity component of Vi os = T4 ys + 
V(R) cos(</>) sin(i) at (i?, </>), where 14 ys is the systemic 
velocity of the galaxy, V{R) is the circular velocity at 
radius R from the galactic center, and <f> is the angle 
from the kinematic major axis of the galaxy. In the sim¬ 
plest case, for a region with a flat rotation curve, i.e. 
V(R) = constant, the observed isovelocity contours of 
the galaxy are given by contours of equal <fi. Therefore, 
velocity gradients exist along any path that perpendicu¬ 
larly crossing those isovelocity contours. The last state¬ 
ment is true for any rotation curve, not just for a flat 
rotation curve. 

On the other hand, we also find velocity gradient in any 
small patches of the data that occupy the GMC regions. 
This velocity gradient can be due to a projection of the 
galaxy rotation, and therefore, can mimic the cloud rota¬ 
tion. In Figure 7, we show the angular momentum vec¬ 
tors of GMCs, overplotted with the isovelocity contours 
of NGC4526. The tendency of the angular momentum 
vectors of the clouds to be tangential to the isovelocity 
contours of the galaxy suggests that the velocity gradi¬ 
ents of the clouds are actually just a projection of the 
galaxy rotation. 

In order to quantify our finding, we create a gas dy¬ 
namical model using t he KinMS (Kinema tic Molecular 
Simulation) package of lDavis et al.l (l2013dh This model 
is basically a purely rotating disk based on a rotation 
curve of the galaxy, i.e. this is what the galaxy look 
like if its dynamic is just a rotation. The rotation curve 
is calcu l ated fr om th e multi-Gaussian expansion (MGE) 
of iCappellaril (|2002f) to the J-band images of HST and 
MDM 1.3-m telescope, and includes the presence of a su - 
permassive black hole at the center ijDavis et al.ll2013bl) . 
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Fig. 7. — The angular momentum vectors of GMCs (black arrows), overplotted with the isovelocity contours of NGC4526 (color coded by 
their projected velocities) convolved with a Gaussian kernel. There is a strong tendency for the vectors to be tangential to the isovelocity 
contours, as expected if the measured velocity gradients of GMCs are just the projection of the galaxy rotation. Correlations between the 
angular momentum of the clouds and isovelocity contours of the galaxy are shown in Figure 8. 




Fig. 8. — Correlations between the model and observed angular speed (left panel) an d the angle of rotation axis 9 (right panel). The 
observed Q and 9 are calculated from the best fit of the velocity field of the cloud (e.g., IGood man et al.M1993l) . The model Q and 9 are 
purely based on the galaxy rotation and created using the KinMS package of IDavis et al.H2013d) as described in the text. The error bars 
are derived from the covariance matrix of the best fit. The excellent one-to-one correlations (solid line) indicate that the velocity gradients 
that we measure are actually just a projection of the galaxy rotation. The dashed lines on the left panel are the standard deviation from 
Monte Carlo simulations that gives the upper limit of the angular speed of the cloud. 
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This MGE fit gives the mass distribution of the galaxy, 
and hence the galaxy rotation curve, parameterized by 
the stellar mass-to-light ratio and the galaxy inclination. 
The model is inclined so that it matches the inclination 
of the observed galaxy (i ~ 79°). Any deviations of 
the data from the model can be caused by small scale 
turbulence, inflow or outflow gas motion, and the cloud 
rotation. 

Then, we measure the angular momentum of the model 
at the location of the observed clouds by using the same 
method as described in §4.1. We find a strong one-to-one 
correlation for both angular speed and rotation angle be¬ 
tween the model and the data, albeit with small scatter 
(Figure 8). The correlation holds true for both SB and 
NSB groups. This reinforces our inference that the veloc¬ 
ity gradients of the clouds are just a consequence of the 
underlying velocity field due to galaxy rotation, i.e. the 
gas within clouds moves following the galaxy rotation. 
Therefore, H that is calculated using the plane-fitting 
method (equation 15) is not the intrinsic angular speed 
of the cloud. 

In this case, the gas within the cloud must rotate due 
to galactic shear. The am ount of shear is given by the 
Oort constant A (|Fleck fc ClarkllTfisTt l: 


Osh ear — 


Av 


V 0 

(dV\ 

A r 

= \2A\ = 

Ro 

Wo 


(16) 


where the subscripts 0 denote the evaluation at the loca¬ 
tion of the GMCs, and V ( R ) is the rotation curve of the 
galaxy. Hereafter, we take this shear as the angular speed 
of the cloud , not fl from equation (15). Furthermore, the 
angle of the rotation axis 9 from the plane-fitting is not 
physically meaningful anymore, because if the cloud ro¬ 
tation is due to galactic shear then their rotation axes 
tend to be parallel to the galaxy rotation axis. 

The intrinsic scatters in Figure 8 may be due to the 
intrinsic angular momentum of the cloud, il c i ou( j, that 
is not originated from the galactic shear. From 1,000 
Monte Carlo simulations, we determine the upper limit 
of the cloud’s angular speed to be H c ioud < 8.6 x 10 -2 
km s _1 pc -1 (dashed line in Figure 8), which is compa- 
rable to the angular speed of GMCs and HI cloud in M33 
dRosolowskv et al1l200.8l llmara et al.ll201 ll) . This upper 
limit is generally smaller than the galactic shear at the 
cloud’s location, i.e. fl c i 0 ud < fl s hea.i, so that H ss fl s hear- 
No correlation between linear and the cloud’s mass is 
found (Figure 8). 

Based on our analysis, any measurement of the ve¬ 
locity gradient of an extragalactic cloud must be per¬ 
formed carefully, to avoid bias due to the projection of 
the galaxy rotation. The only exception is if the galaxy 
is nearly face-on (i m 0°), as the line of sight velocity due 
to galaxy rotation is negligible. 


4.3. Stability of Rotating Clouds 

In §3.2, we assumed that the clouds were supported 
by turbulence only. In the presence of rotation (due 
to shear), the rotational energy also contributes to the 
clouds’ stability against gravitational collapse. Here, we 
define the parameter 7 rot as the ratio of rotational over 
turbulent energy: 


Trot — 


pfl 2 R 2 

3a 2 


(17) 



log Trot 


Fig. 9. — Distribution of 7rot, defined as the ratio between ro¬ 
tational and turbulent energy. The dashed line is the mean of the 
log-normal fit (red). For the majority of the clouds, the rotational 
energy due to galactic shear is smaller than the turbulent energy. 
From the correlation between 7 ro t and luminous mass (shown as 
an inset), we infer that the relative importance of rotation over 
turbulence is increasing for more massive clouds. 



Fig. 10. — Same as Figure 9, but for /3 ro t, which is defined as 
the ratio between rotational and gravitational energy. Most of the 
clouds are not in rotational equilibrium, i.e. /3 ro t < 1. The four 
clouds in the tail of the distribution are the innermost clouds that 
suffer strong shear. From the correlation between /3 ro t and lumi¬ 
nous mass (shown as an inset), we infer that the relative importance 
of rotation over gravity is increasing for more massive clouds. 


where y = 2/5, th e value for a uniform sphere 

(Goodman et al.ill99.f b 

In Figure 9, we show the histogram of the 7 ro t values of 
the clouds. A log-normal fit to the distribution yields a 
mean 7 rot ~ 0.24 with a standard deviation of 0.37 dex. 
Roughly 92% of the resolved clouds have 7 ro t < 1. This 
means that the rotational energy is smaller than the tur¬ 
bulent energy for the vast majority of the clouds. There 
is also a correlation with mass as 7 ro t oc M 1 ^ 3=t0 ' 16 , so 
that the relative importance of rotation over turbulence 
is increasing for more massive clouds. 

It is also useful to define the ratio between rotational 
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Fig. 11.— Same as Figure 9, but for the virial parameter a. 
The dashed line is the mean Gaussian fit, and the dotted line is 
a = 1.13. The four clouds in the tail of the distribution are the 
innermost clouds that undergo a strong shear. Except for those 
four clouds, the cloud population is in gravitational equilibrium, 
i.e. the mean a ^ 1.26. The inset shows anti-correlation 
the virial parameter and the luminous mass. 

kinetic energy and self-gravitational energy: 

_ 1 p 0 2 P 3 
^ rot “ 2 q GM ’ 

where p/ g = 2/3, the value for a uniform 
(Goodman et al.l 19931 ). Clouds with /3 rot ft 
rotationally-stable against gravitational collapse. 

96% percent of the resolved clouds have /? < 1. 
normal fit to the distribution yields a mean /? rot 
with a standard deviation of 0.36 dex (Figure 10). This 
means that the rotational energy is smaller than the 
gravitational energy for the vast majority of the clouds. 
Furthermore, there is a correlation with mass as /3 ro t oc 
-^ium° ±0 21 > so that the relative importance of rotation 
over gravity is increasing for more massive clouds. 

Finally, the virial parameter a, which includes turbu¬ 
lence, gravity, and rotation, can be expressed as 

ot = /3 r ot (1 H-^ • (19) 

V 7rot J 

Non-magnetic, rotating clouds with a ft 1.13 are in 
virial-equilibrium (BM92), while clouds with a>l are 
either pressure-confined clouds, or still gravitationally 
bound but with an underestimated CO-to-H 2 conversion 
factor. 

The distribution of a values is shown in Figure 11, 
where it can be approximated by a log-normal distri¬ 
bution with a tail at the high end. The log-normal 
mean is a ft 1.26 with a standard deviation of 0.15 dex. 
This suggests that the GMC population in NGC4526 is 
gravitationally-bound, even after the inclusion of rota¬ 
tional energy. If we define non-gravitational equilibrium 
as having a > 3.5 ft 0.5 dex (i.e. the tail of the dis¬ 
tribution), then only ft! 4% of resolved clouds are not 
gravitationally-bound. As shown in the inset in Figure 
11, by excluding those four clouds, a has a shallow cor¬ 
relation with mass as a oc A/ 1 “^ 19±0 H5 . 
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Fig. 12.— Top: ratio of internal and gravitational pressure as a 
function of distance from the galactic center. Except for the four in¬ 
nermost clouds, the data points are consistent with gravitationally- 
bound clouds. The lcr scatters are indicated with dashed lines. The 
SMBH radius of influence (ft 45 pc) are indicated as vertical dotted 
lines. Bottom: galaxy circular velocity curve. The sharp increase 
near t h e center is d ue to the presence of a supermassive black hole 
(Davis et alJl2013b ). 

Further investigation reveals that the four 
gravitationally-unbound clouds (with log a > 0.5) 
are the clouds closest to the center of the galaxy, at 
a distance of ft 10, 34, 42, and 54 pc. These clouds 
suffer strong shear (|Au/Ar| ~ 10 km s -1 pc -1 ) due to 
the presence of a supermassive black hole (SMBH) at 
the center of NGC4526, with Mrh ft 4.5 x 10 8 M ^ and 
radius of influence of ft 45 pc dDavis et al.ll2013bT) . This 
SMBH makes the circular velocity curve of the galaxy 
increases abruptly towards the galactic center (Figure 
12), and hence yields a large Oort A constant. 

5. DISCUSSION 
5.1. Pressure Balance 

In general, the dynamical equilibrium state of a cloud 
can be written as 

-Pint + Pb = Pg + Pextj (20) 

where Pj n t ft Ptm-b (1 + Trot) is the internal pressure of 
the cloud, including the correction factor (1 + 7 rot ), due 
to the contribution of rotation. Pturb = per 2 is the kinetic 
pressure due to turbulence, P ex t is the external pressure 
of the ambient medium, and Pb = B 2 /8t: is the magnetic 
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pressure. Here, we assume that the thermal pressure is 
much smaller than the pressure due to turbulent motion, 
and it is neglected. Pg is the internal gas pressure that 
is required to support the cloud against gravity in the 
absence of any other forces (BM92): 


Pg 

k 


1.3 <Fg 




( 21 ) 


where < j>G is a dimensionless factor that measures the 
ratio between the gravitational pressure of ellipsoidal and 
spherical clouds, and depends only on the cloud’s axis 
ratio: cr ma j/cr m i n for prolate clouds and cr m i n /cr ma j for 
oblate clouds. Here, we assume all clouds are prolate. 
The value of < 4>g > for oblate clouds is within the 
uncertainty of that for prolate clouds. 

For non-magnetic, rotating, gravitationally-bound 
clouds: Pint ss Pg- In Figure 12, we plot log(P int /P G ) 
versus the distance of clouds from the galactic center. 
Except for the four innermost clouds, log(P; n t/PG) has a 
mean value of « —0.03 and standard deviation of ss 0.18 
dex, consistent with a gravitationally bound state. We 
do not include the error bars of individual data point to 
calculate those values. 

As we mentioned in §4.3, there are four central clouds 
that experience strong galactic shear. By using equation 
(20) and an assumption of zero magnetic pressure, the ex¬ 
ternal pressure that is required to bind the clouds against 
galactic shear is P ex t ~ 10 9 K cm -3 , which is extremely 
high and unrealistic. We argue that these are unbound 
clouds. If nothing balances the shear, then these clouds 
will be ripped apart by strong shear within a timescale of 
~ 2-7 t/| Au/Ar| ~ 1 Myr. This timescale is smalle r tha n 
the e xpected lifetimes of GMCs (rs 30 Myr; iBlitz fc Shul 
11980) based on the dumpiness nature of GMCs and de¬ 
struction processes from massive star formation inside 
GMC complexes. Another evidence of the role of shear 
in the destruction of GMCs is found in the M51 disk, 
where fMivamoto et al. 1 2014 ) reported that the locations 
of giant molecular associations are anti-correlated with 
the shear strength. 


5.2. Does a Size-linewidth Relation Really Exist? 

Clouds in NGC4526 do not show a size-linewidth rela¬ 
tion (see Figure 5), in contrast to the previous arguments 
that supported the Larson’s ‘law 1 (e .g. in the Milky Wa y 
disk; S87, and Local Group galaxies: lBolatto et ahfeOQS) . 
Most of the Local Group members are late-type galaxies, 
while our clouds reside in the bulge of a SO-type galaxy. 
Does this discrepancy suggest that GMC properties in 
the early-type galaxies are intrinsically different? Here, 
we argue that the size-linewidth relation may not exist 
in all galaxy morphologies, because the cloud’s size and 
linewidth within a single galaxy only have weak to mod¬ 
est correlation coefficients. 

Recent studies of GMCs in spirals, such as M3 3 
dGratier et all 12012) and M51 ([Colombo et all 12014) . 
also found no clear size-linewidth relation, with a Spear¬ 
man correlation coefficient r sp of 0.12 and 0.16, respec¬ 
tively. A mo dest c orrelation (r sp = 0.51) was found by 
iHever et al.1 (12009) . who re-examined S87 clouds in the 
Milky Way disk using more sensitive instruments. Fur¬ 
thermore, clouds in the LMC also show a weak correla - 
tion (r sp = 0.37: iHuehes et all 120131 : iWong et all 12011) . 
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Fig. 13.— Plots of size vs. velocity dispersion for extragalac- 
tic GMC populations. The contours enclose 68% of the distribu¬ 
tion of data points of a given galaxy: NGC4 526 (black; this pa¬ 
per), Milky Way disk (blue; IHever et aT1l2009f ). Milky Way center 

S ; [Oka ct al. 1998), Large Magellanic Cloud (red; IWong et abl 
, M51 (green ; IColombo et al.l 1201411 . and M33 (magenta; 
er et all 12012 ;). The centers of a Gaussian fit to each dis¬ 
tribution are shown as filled circles. The dashed line is the Milky 
Way disk relatio nshi p (S87 ), the dotted line is the Milky Way cen¬ 
ter relationship ( Oka et al.l 119981. and the solid line is the Local 
Group relationship iiBoIatto et al.l[2008 ). 


These evidence suggest that the size-linewidth relation 
may not exist in all galaxy morphologies. 

If the argument for a size-linewidth relation is not con¬ 
clusive for GMCs within a single galaxy, then how about 
a comp i lation of JdMC data from various galaxies (e.g. 
iBolatto et al.ll2008) ? In this case, one must pay attention 
to different data sets that have different physical resolu¬ 
tions and sensitivities. Coarse resolution and low S/N 
observations can only measure average properties within 
a larger area, without the ability to decompose the CO 
structure into multiple smaller clouds, while finer reso¬ 
lution observations tend to over-decompose CO emission 
into smaller scale structures. This means that the iden¬ 
tified GMCs in different data sets are likely to probe 
different scales of CO emission. This bias, which is due 
to the ability to decompose structure in GMCs, is sep¬ 
arate from the bias of measured properties due to finite 
resolution and sensitivity, which has been minimized by 
the CPROPS program. Hence, plotting those data in 
a size-linewidth diagram leads us to co m pare d i fferent 
structure of GMCs. For example, IBolatto et al.l ()2008) 
compare composite extragalactic GMCs that have been 
observed with a range of resolution from ~ 6 pc (about 
the size of a clump) to ~ 117 pc (about the size of a gi¬ 
ant molecular association). Therefore, any scatter in the 
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size vs. line width diagram is overcome by the large range 
of GMC size, which gives rise to a slope. Even in this 
case, however, the correlation coefficient is still moderate 
(r sp =0.57). 

A self-consistent study of the size-linewidth relation, 
then, requires a common physical resolution and sen¬ 
sitivity across the extragalactic GMC data sets. Any 
similarity or discrepancy among the extragalactic GMCs 
measured in th is wa y wo uld then be genuine. In a recent 
work, iHughes et al.l ( 2013 1 showed that a size-linewidth 
relation is apparent when M51, M33, and LMC data were 
analyzed at the ir original (different) resolutions and sen¬ 
sitivities (as in iBolatto et al.l 120081) , but no compelling 
evidence was found when the data were degraded to a sin¬ 
gle (‘matched 1 ) common resolution and sensitivity. From 
the ‘matched 1 data, they infered that GMCs in M51 are 
in general larger, brighter, and have higher velocity dis¬ 
persions than equivalent structures in M33 and the LMC, 
which can be interpreted as a genuine variation of GMC 
properties. 

In t his respect, we can compare NGC4526 and LMC 
data (jWong et al.l l2011f ). since they have comparable 
physical resolutions and sensitivities (~ 20 pc and cr rms ~ 
0.7 K for NGC4526, and ~ 11 pc and cr rms « 0.3 K 
for LMC). We find that GMCs in NGC4526 tend to 
have higher velocity dispersions than equal size clouds 
in the LMC (Figure 13). With respect to the Milky Way 
disk, clouds in NGC4526 lie above the Milky Way disk 
size-linewidth relation (S87) by a factor of ~ 3 (Fig¬ 
ure 13) and above the Milky Way size-luminosity rela¬ 
tion by a factor of ~ 5 (Figure 5), which implies that 
clouds in NGC4526 are brighter and more turbulent than 
similarly-sized clouds in the Milky Way. In contrast, 
NGC4526 clouds are less turbulent than the Galactic 
center clouds by a factor of ~ 0.4 dex (Figure 13). This 
genuine variation of GMC properties may be influe nced 
by different environment between galaxies (IHughes et al.1 
I2013D and is discussed in §5.4. 

5.3. Variations of GMC Surface Density 

If the standard size-linewidth relation ( cr v oc A 1 / 2 ) is 
valid for GMCs, then as a consequence the mass of grav¬ 
itationally bound clouds is M v j r oc R 2 (c.f. equation 4), 
and hence the mass surface dens i ty E = M^/ttR 2 = 
constant. However, iHever et al.l (l2009t) . who revisited 
the GMCs of S87 using more sensitive and better sam¬ 
pled data, found that the surface density is actually not 
constant, and the coefficient of the size-linewidth rela¬ 
tion (Co = cr„i? -1 / 2 ) correlates with the surface density 
as Co oc E 1 / 2 . This relation is expected from gravita¬ 
tional equilibrium (equation 4), and does not depend 
on whether the clouds follow the size-linewidth relation 
or not. The same relation (Co oc E 1 / 2 ) also holds for 

pressure equilibrium, but with E = E c oc P e f t (Field 
et al. 2011), where E c is the critical surface densit y 
of a Bonnor-Ebert sphere dBonnod 119561 : lEbertl H9551) . 
The difference between the two is that pressure equilib¬ 
rium has a higher normali zatio n than the gravitational 
equilibrium. The IHever et al.l (2009) data favor pres- 
sure e quilibrium rather than gravitational equilibrium 
(IField et alJl20lll) . 

In Figure 14, we compile extragalactic GMC data. The 
contours enclose 68% of the distribution of data points of 



Fig. 14.— Correlation between a v R~^C anc [ surface density 
for extragalactic GMC populations (color coded as in Figure 13). 
The contours enclose 68% of the distribution of data points of a 
given galaxy. The solid ‘V‘-curves are the pressure equilibrium 
condition of an isothermal sphere for various external pressures 
(P ex t/k = lO 7 ,10 6 ,10 5 ,10 4 , and 0 K cm -3 , respectively, from top 
to bottom; IField et alJI20lflD . The dashed line is the locus of crit¬ 
ical surface density for a Bonnor-Ebert sphere. The dotted line 
is the gravitational equilibrium of a constant density sphere. This 
plot shows that the surface density of GMCs is not constant as pre¬ 
viously believed. Inset: the surface density function of the clouds 
in NGC4526. 

each galaxy. Our compilation of extragalactic GMC data 
shows that the surface density is not constant, but varies 
from ~ 10 to 3000 M@ pc -2 . GMCs in the Milky Way 
disk, LMC, and M33 have lower surface densities than 
GMCs in M51 and NGC4526. The median surface den¬ 
sity of the NGC4526 clouds is E me d ~ 1.2 x 10 3 M© pc -2 , 
which is ~ 7 times greater than in the Milky Way disk 
clouds (170 M© pc -2 ; S87). However, NGC4526 clouds 
have similar surface density as the Galactic center clouds. 
This high surface density may be a common feature for 
clouds in the galaxy bulge. 

In Figure 14, we also see that there is a correlation 
between cr V R~ 1 / 2 and surface density, as expected from 
gravitational (and pressure) equilibrium. The Milky Way 
disk, LMC, M33, M51, and NGC4526 clouds roughly 
follow a v R~ x ! 2 oc E 1 / 2 , but the Milky Way, LMC, 
and M33 clouds have higher normalizations, i.e. they 
lie above the gravita tional equilibrium relation (dotted 
line). As IField et al.l (|2011l) suggested, Milky Way disk 
clouds are likely to be in pressure-equilibrium (dashed 
line) rather than gravitational equilibrium, and hence 
have a higher normalization factor. With the exception 
of Galactic center clouds, it is interesting to note the 
trend that GMCs with lower surface densities tend to 
be in pressure equilibrium. The Galactic center clouds 
are unique because they are pressure-bound clouds with 
M v i r ~ 10 Mi um and they reside in th e hi gh external pres¬ 
sure environment (Ok a et all 119981 iMivazaki fe Tsuboil 
( 200 ( 1 ) . 

Based on this finding (Figure 14), we argue against 
the current mainstream view regarding the constancy of 
cloud surface density. Indeed, theoretical studie s (e.g. 
lKegellll989t iBallesteros-Paredes fc Mac Lowll2002D found 
that limited observational sensitivities can give biased re- 
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suits, so that previous measurements of the surface den¬ 
sity of GMCs appear constant. Moreover, size, veloc¬ 
ity dispersion, and surface density are correlated with 
each other as expected from gravitational equilibrium 
(or pressure equilibrium as in the Milky Way, LMC, and 
M33). 

Here, we propose a modified version of the Larson’s re¬ 
lations. (1) The clouds are in either gravitational or pres¬ 
sure equilibrium; the relative contribution of gravity and 
external pressure to cloud stability need further study 
on a cloud-by-cloud basis. (2) The coefficient of the size- 
linewidth relation depends on the cloud surface density 


as Co oc E 1 / 2 . This relation also holds true in the pres¬ 
sure equilibrium case, where surface density depends on 

1/2 

the extern al pres sure of the ambient medium (E oc P e / t ; 
[Field et alJl20H . (3) The cloud surface densities are 
not all the same, but may depend_on_the environment, 
such as the external pressure ( Elmegreeiil 119931) , inter¬ 
stellar radiation field strength ( McKee[ 1989fl. and inter- 
stellar gas flow and turbulence ( Vazauez-Semadeni et al.1 
I2007T) . which need further investigation. These environ¬ 
mental dependencies could explain variations o f GMC 
properties across differen t galaxies (I Hughes et al.l I 20 T 0 I : 
IRosolowskv fc Blitzl 120051) . 


5.4. Environmental Effects 

GMCs in NGC4526 are denser and more turbulent 
than those in the Galactic disk, but have similar sur¬ 
face density and less turbulent than those in the Galac¬ 
tic center (Figure 13 and 14). These differences may be 
caused by different environment, such as the interstellar 
radiation field strength (ISRF) and the external ambient 
pressure (P ex t), between NGC4526 and the Milky Way. 
Here, we infer the ISRF and P ex t based on the global 
properties of the galaxy and discuss their possible roles 
to explain the differences between GMCs in NGC4526 
and the Milky Way. 

5.4.1. Comparison with GMCs in the Milky Way disk 

ICiesla et al.l (|2014T) have used the Herscliel photometric 
data to derive the dust spectral energy distribution of 322 
nearby galaxies, including NGC 4526. Th ey f it the data 
with the dust emission model of IDraine fe Lil (120071) . In 
this model, a large fraction of dust is located in the diffuse 
interstellar medium, exposed to a single ISRF with inten¬ 
sity U = U min • We define U as the intensity normalized 
to the Milky Way value, i.e. ISRF = U x IS RFmw, where 
ISRF mw is the ISRF of the Milky Way (iMathis et al.1 
119831) . In addition, there is a small fraction ( 7 ) of dust 
located in regions where the ISRF is more intense (e.g. 
photodissociation regions), with ISRF rangin g from I/m in 
to U m ax and described by a power law U~ a . iDraine &; Lil 
(l2007ll found t/ max = 10 6 and a = 2 are the best fit 
to the SINGS sample (|Kennicutt et al.l[20031) . The free 
par ameters of the model are then reduced to t/ m ; n and 
70 ICiesla et al.1 (12014 1 found the best fit parameters of 
U min = 3.92±0.32 and 7 = (0.19±0.09)% for NGC_4526. 

Then, we can calculate the mean ISRF (U) in 
NGC4526, weighted by the dust mass, by using equa- 

7 There is a third parameter of the model, namely the fraction 
of dust mass contributed by Polycyclic Aromatic Hydrocarbons 
(PAHs), but we do not need it to calculate the mean ISRF. 


tion (17) of IDraine fc Lil (120071k 


U = (1 - 7) Pmin ± 7 « 4.0. ( 22 ) 

^min c/ma X 

Thus, the ISRF in NGC4526 is ~ 4 times higher than in 
the Milky Way. 

This higher ISRF means higher photoionization rate 
to the molecular gas, so that the CO emission emerges 
from a deeper layer within the cold gas, i.e. at higher ex¬ 
tinction (Ay) than in the typical Milky Way disk. This 
higher extin ctio n tran slates into higher gas surface den¬ 
sity via (lMcKeelll989l) 

Eh 2 = 22.3 ^ Mq pc" 2 , (23) 


where <5 gr is the ratio of the extinction per hydrogen 
nucleus in the cloud to the standard value given by 
ISnitzed (|1978l k Generally, <5 gr = 1 in the Milky Way. 
We do not know the value of S gI in NGC4526, but we 
can assume th at is proporti onal to the metallicity as 
<5 gr ~ Z/Zq (|Bolatto et al. 120081). By using equation (23) 
and log (Z/Zq) ss 0.2 ( Davis et al.ll2013a)k we can esti¬ 
mate the value of Ay that is required to reproduce the 
observed surface density of GMCs in NGC4526. Cloud’s 
surface density in NGC4526 is ~ 7 times higher than 
those values in the Milky Way disk (Figure 14), so that 
Ay in NGC4526 is estimated to be 4.4 times higher than 
that in the Mil ky Way disk. 

Furthermore, IMcKeel (|1989l ') predicts that the velocity 
dispersion of GMCs is proportional to the square-root of 


Ay and the cloud size as a v oc (Ay/6 gr) 1 / 2 R 1 ^ 2 . This 
relation arises naturally from the gravitational equilib¬ 
rium state and by using Ay as a proxy of surface density 
as in equation (23). By using the estimated value of Ay 
above, the theory predicts that the velocity dispersion 
in NGC4526 clouds is about 2.1 times higher than the 
velocity dispersion of equal-size clouds in the Milky Way 
disk. This prediction is in agreement with our measure¬ 
ments, which show the velocity dispersion of equal-size 
clouds in NGC4526 is higher than those in the Milky 
Way disk by a factor of ~ 3 (Figure 5 and 13). Thus, we 
speculate that the surface density and velocity dispersion 
in NGC4526 clouds are higher because those clouds have 
higher extinction than clouds in the Milky Way disk. 


5.4.2. Comparison with GMCs in the Galactic center 

Clouds in NGC4526 have similar surface density and a 
smaller velocity dispersio n than equal size clouds in the 
Galactic center (|Oka et al.l 119981) by a factor of ~ 0.4 
dex (Figure 13 and 14 ). This may be due to the fact 
th at jOka et al.l (Il998f) clouds and iMivazaki fc Tsuboil 
(I200(fr clumps are in pressure equilibrium, rather than 
gravitational equilibrium as in NGC4526 clouds. From 
the equation (20), by neglecting the magnetic pressure 
term, this means the internal pressure of pressure-bound 
clouds needs to balance against gravity and external pres¬ 
sure, while gravitationally bound clouds needs to balance 
against gravity only. Therefore, for a given cloud mass 
and radius, the velocity dispersion of pressure-bound 
clouds is higher than the velocity dispersion of gravita¬ 
tionally bound clouds in order to maintain a dynamical 
equilibrium state. 
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We do the following calculations to support our argu¬ 
ment. By neglecting the magnetic field, gravitationally 
bound clouds have a 2 vir = Pg p _1 , while pressure-bound 
clouds have a 2 s = ( Pg + P ex t) p _1 , where p oc £ i? _1 
and Pq oc £ 2 (BM92). For equal-size clouds with similar 
density, the ratio between the two is 

/ p \ 1/2 

^pres = / 1 + £«ct \ ( 2 4) 

&v,v ir V Pg ) 


For spherical clouds, like those in NGC4526, with mass 
~ 10 6 Mq and radius ~ 20 pc, Pg is ~ 8 x 10 6 K cm -3 . 
The external pressure P ext hr the Galactic center is rather 
uncertain. For P ext /fc between ~ 5x 10 6 ( Suergel fc Blit3 
I1992T) and ~ 1 x 10 s (|Mivazaki fe Tsuboil l2000h . we get 
the Galactic center clouds are expected to have higher 
velocity dispersion by a factor of 1.3 to 3.7 (c.f. equation 
24). This range is also in agreement with our measure¬ 
ment (~ 0.4 dex « 2.5; Figure 13). 

If this is true, then why are the clouds in NGC4526 
bulge in gravitational equilibrium but clouds in the 
Galactic center pressure bound? To get insight into 
this question, we estimate the global ambient hydrostatic 
press ure Ph in N GC45 26 as a proxy of the external pres¬ 
sure (lElmegr eenl 1 19891 1: 


Ph = — (llg + ^ £* 

2 V er* 


(25) 


where £ g = Eri is the ambient gas surface density, £* is 
the stellar surface density, a g is the ambient gas velocity 
dispersion, and cr* is t he stellar velocity dispe rsion. 

We use ATLAS-3D dCappellari et all 1201 111 results to 
get the stellar properties of NGC4526: cr* « 233.3 
km/s (at the ce ntral 1 kpc), L m 3.13 x lO lo P 0 . r 
( Caggellari^^L 2013aJ), stellar M/L r « 5.6 M 0 /L 0 ,r 
( Cappellari et al. 2013blb and R e ~ 74.1” w 5.9 kpc 
( Krainovic et al.1 120131b The quoted luminosity value 
is global, so the enclosed luminosity at R e is simply 
half the quoted value, i.e. L(R e ) ss 1.57 x lO 10 L 0 . 
Thus, the stellar surface density can be estimated as 
£* « M(R e )/nR 2 e « 804 M & pc" 2 . 

HI is undetected in NGC4526 with uppe r limit of 
Mm < 1.9 x 10 7 Mq (|Lucero fc Yound [20131. The lin¬ 


ear resolution of their observations is 5.4 x 4.2 kpc 2 . 
This gives an upper limit of the gas surface density as 
£ s < 0.27 Mq pc -2 . We do not know the value of a g , 
so we assume a g = 10 km/s (jBlitz fc Rosolowskvll2006l h 
Taken all together, we estimate the external pressure 
to be Ph/k < 1.5x 10 6 K cm -3 . Thus, unlike the Galactic 
center clouds, where P ex t ~ Pg, NGC4526 clouds have 
Pext < Pg- This may cause the Galactic center clouds 
are pressure bound, while clouds in NGC4526 are grav¬ 
itationally bound (since the external pressure is small 
with respect to the gravity). This small external pres¬ 
sure is presumably due to lack of HI in NGC4526, which 
may be caused by ram pressure or hot gas evaporation 
as the galaxy resides in the Virgo cluster. 


6. SUMMARY 

We identify 241 GMCs based on 12 CO(2 — 1) observa¬ 
tion at ~ 20 pc resolution in the galaxy NGC4526 using 
the CPROPS program (RL06), where 103 of them are 
spatially resolved. As a population, the clouds are in 


gravitational equilibrium. A log-normal fit to the pop¬ 
ulation yields a mean virial parameter a ~ 0.99 with a 
standard deviation of ~ 0.14 dex. 

The cloud mass distribution follow dN/dM oc 
M“ 2 - 39± o o 3 ; steeper than in the inner Milky Way but 
comparable to what others find in several other galaxies. 
Since the exponent is less than —2, the total molecu¬ 
lar mass is dominated by the contribution of low mass 
clouds. The data favor a truncated-distribution with 
truncation mass of 4.12 x 1O 6 M 0 . 

In general, clouds in NGC4526 are more luminous and 
more turbulent than equal-size clouds in the Milky Way 
disk by a factor of ~ 5 and ~ 3, respectively. Moreover, 
the surface density of GMCs in NGC4526 is ~ 7 times 
higher than those in the Milky Way disk. These differ¬ 
ences may be due to higher ISRF and cloud extinction 
{Ay) in NGC4526, so that the CO emission emerges from 
a deeper layer in the cold gas, and hence, a higher gas 
density region. 

On the other hand, NGC4526 clouds are less turbulent 
than the Galactic center clouds. This may be caused 
by different equilibrium state of GMCs: Galactic cen¬ 
ter clouds are pressure-bound, while clouds in NGC4526 
are gravitationally bound. The velocity dispersion of the 
pressure-bound clouds needs to balance both gravity and 
the external pressure, while gravitationally bound clouds 
just need to balance gravity only. Indeed, our estimation 
shows that the external pressure in NGC4526 is smaller 
than the gravity, so that the external pressure is less im¬ 
portant in the dynamical state of NGC4526 clouds. This 
situation is different in the Galactic center, where the 
external pressure is comparable or higher than the cloud 
self-gravity. 

We find no size-linewidth correlation in NGC4526 in 
contrast to what is expected from Larson’s relation. This 
finding is robust against the choice of the input param¬ 
eters of the CPROPS program or different measurement 
methods (Appendix D). This implies that the surface 
density of GMCs is not constant, but follows the relation 
avR- 1 ' 2 oc E 1 / 2 as expected from gravitational equilib¬ 
rium. 

In the kinematic analysis, we find that the velocity gra¬ 
dient of individual clouds are just a consequence of galac¬ 
tic rotation. In this case, if the clouds are rotating, then 
the rotation follows the galactic shear given by the Oort 
A constant at the location of the cloud. We calculate 
9 s hear and find that 92% of resolved clouds have a tur¬ 
bulent energy exceeding the rotational energy, and 96% 
of resolved clouds have a gravitational energy exceeding 
the rotational energy. This means rotational energy is a 
minor contribution to the clouds’ dynamical stability. 

Even with the inclusion of rotational energy, the cloud 
population is still in gravitational equilibrium. The dis¬ 
tribution of the virial parameter can be approximated 
by a log-normal distribution with a tail at the high-end. 
The mean of the distribution is a ~ 1.26 with a stan¬ 
dard deviation of ~ 0.15 dex. There are only 4 clouds 
with a > 3.5. These clouds are the innermost clouds and 
undergo extreme galactic shear. These unbound clouds 
should be ripped apart in a timescale of less than ~ 1 
Myr. 
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APPENDIX 


A. CLOUD IDENTIFICATION ALGORITHM 

We use the modified CLUMPFIND algorithm dWilliams et al.1 119941) . implemented in the CPROPS program 
(Rosolowskv & Leroy 2006, hereafter RL06) to identify GMCs in NGC4526. Below are descriptions of the code, 
together with our chosen values of input parameters. 

First, the program identifies connected regions of significant emission as islands. An island is defined as CO emission 
that has at least one pixel higher than 3cr rms (the threshold parameter of the program) and extends to all connected 
pixels with emission higher than 2er rms (the edge parameter of the program). An island consists of one or more clouds 
after the decomposition process. We set the minimum volume of islands to be 1 beamwidth 2 x 1 velocity channel. We 
choose these values to include any possible small island in our data, since our resolution is somewhat comparable to 
the typical size of Milky Way GMCs. 

The decomposition of each island begins by looking for local maxima. Local maxima are identified by looking for 
pixels that are greater than or equal to all neighbors in a 1 beamwidth 2 x 1 velocity channel volume. Our choice of 
these values is to separate an island into potentially smaller clouds. 

For each local maximum, working from the local maximum that has the lowest emission to the highest, the algorithm 
identifies pixels associated exclusively with each local maximum by contouring the data cube in three dimensions. If 
the emission of a local maximum is less than na rms above the merge level with neighboring local maxima, or there 
are fewer than m x beamwidth pixels associated with the local maxima, then the local maximum is removed from 
consideration. The merge level is the contour level that encloses two neighboring local maxima. The purpose of this 
decimation process is to remove spurious peaks of noise, i.e. false clouds. Higher values of n and m gives a smaller 
probability of false clouds, at the cost of losing small genuine clouds. We adopt m = 0.5 to account for small clouds 
and determine the best value of n from simulations (described in Appendix B). 

Then, the program decides whether two neighboring clouds are a merged cloud or distinct clouds. The algorithm 
compares the values of emission moments for the separated and combined clouds. If the flux F and moment er of an 
individual cloud differ by more than a fraction of the flux and the moment of the merged cloud ( SF/F and Sa/a), 
then the local maxima are categorized as distinct. The values of 5F/F and Sa/a are chosen from the simulations 
(Appendix B). 

At the end of the process, the program calculates the properties of the clouds (described in §3.1) and records them 
into a catalog (Table 1). 


B. DECOMPOSITION PARAMETERS 

In order to choose the best decomposition parameters (n, SF/F , and Sa/a) that are suitable for our data, we create 
two simulated clouds with a 3D-Gaussian shape in a single data cube, and add the typical noise of our observations. 
The dispersions of the Gaussian are a x = a y = 1 beam width and a v = 1 velocity channel, so that they are resolved by 
the antenna beam and spectral channels. We vary the separation (center-to-center) of those gaussians in units of R, 
which is defined as R = 2er x . Two Gaussians are almost fully resolved if they are separated by a distance larger than 
2 R. We also vary the peak S/N of the Gaussians from 3er rms to 6er rms , to take into account any possible dependence 
of our simulations on S/N. We run 10 simulations for each choice of separation distance and peak S/N, so that the 
results are statistically robust. 

We feed the simulated data cubes into the CPROPS program and vary the three input parameters that drive the 
decomposition of the clouds (n, SF/F, and Sa/a) from 0 to 3 with an increment of 0.5. The program then identifies 
the number of clouds in a given data cube. For various decomposition parameters and peak S/N, we plot the average 
number of clouds identified by the program against the separation distance in Figure 15. The program successfully 
resolves two clouds for a separation distance larger than 2 R. However, for blended Gaussians (separated by a distance 
shorter than 2 R), the values of n = SF/F = Sa/a = 1 best recover the correct number of clouds at all S/N. Therefore, 
we adopt these values as the decomposition parameters for our data. In Appendix D, we further show that the results 
of our studies are not sensitive to the choice of decomposition parameters. 
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Separation [R] 


Fig. 15. — Average number of recovered clouds as a function of the separation (center-to-center) between two clouds. The separation 
distance is in units of R = 2<r x , where cr x is the dispersion of the Gaussian clouds. Different panels show different peak S/N of the simulated 
clouds. In each panel, different values of the decomposition parameters (n, 5F/F , and Safer) are shown in different colors. The correct 
number of clouds is indicated by the dashed lines. The program recovers approximately two clouds for separations of at least 2 R. Unity 
input parameters approach the correct number of clouds even for blended clouds. 


C. PROBABILITY ANALYSIS OF REAL DETECTIONS 

In o rder to chec k t he p robability that the identified clouds are real, we do a probability analysis similar to that given 
in lEneargiola et ahl (120031) . If we have n-adjacent channels with the same brightness temperature T [, = ka rms , then the 
probability of this being a false detection is P n {k) = [0.5 x erfc(fc/-\/2)] n , where erfc is the Gaussian complementary 
error function. The probability of real detection is P rea i = 1 — N tr i a \Pn(k) for N tr i a \Pn(k ) <C 1. Here, Aerial = A'/n, 
where N ss 2.37 x 10' is the number of pixels in our data cube. If the pixels are not independent due to beam 
convolution and spectral smoothing, then the inferred P rea i is smaller than the true P rea i- We set the edge parameter 
of the CPROPS program to be 2<r rms , so that all pixels in a cloud must have Tj, > 2er rms . Hence, the probability that 
a cloud occupying n-adjacent channels is a real detection is P rea i > 1 — A tr i a iP rl (2cr rms ). In this case, P rea i > 0.97 for 
n = 5. The smallest identified cloud has total number of pixels n = 13. This suggests that it is unlikely that we detect 
false clouds, so we assume all identified GMCs are real structures. 

D. CHECKING BIAS AGAINST THE CHOICE OF INPUT PARAMETERS 

We evaluate how our results are affected by the choice of input parameters of the CPROPS program (as described 
in Appendix B) and the methods to measure the cloud properties (as described in §3). We rerun the program and 
recalculate the cloud properties (radius, velocity dispersion, and luminosity) using various decomposition parame¬ 
ters ( n,6F/F , and ^cr/cr), edge and threshold parameters for islands, and methods of measurement (by excluding 
deconvolution, excluding extrapolation, and excluding both deconvolution and extrapolation). 

The decomposition parameters are varied from 1 to 3 with unity increment (first column of Figure 16), the edge 
parameter varies from 2<j rms to 3<r rms , and the threshold parameter varies from 3<r rms to 4cr rms (second column of 
Figure 16). Different methods of measurements are given in the third column of Figure 16. Each row in Figure 16 
shows the plots between various cloud properties, i.e. size vs. linewidth in the first row, linewidth vs. luminosity in 
the second row, and size vs. luminosity in the third row. Larson’s relations for the Milky Way are shown as the solid 
lines. For comparison, the dashed lines are Larson’s relations with different normalization factors: 3, 0.05, and 5 for 
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Fig. 16.— Plots of GMC properties (radius, velocity dispersion, and luminosity) for various decomposition parameters (left), islands 
parameters (middle), and methods of measurement (right). The black lines are Larson’s relations and the dashed lines are Larson’s 
relations with higher or lower normalization factors. There is no obvious bias due to the choice of input parameters of the program. 


R vs. a v , a v vs. L , and R vs. L , respectively. To check how strong the correlation between various cloud properties 
are, we tabulate the Spearman correlation coefficients r sp in Table 3. 

For each plot, we build an estimate of the probability density function (PDF) based on the data scatter in two- 
dimension, using the kernel-density-estimate method in Scipy. The contours in Figure 16 enclose to the 68% confidence 
level of the PDFs. We also fit the PDFs with a 2D-Gaussian and show the Gaussian centers as filled circles in Figure 
16. The center and dispersion of the Gaussians are tabulated in Table 4. For clarity, we do not show the data points. 
In addition, we check whether the PDFs in a given panel of Figure 16 are sampled from the same parent distribution 
or not. This is tested with the two-dimensional Kolmogorov-Smirnov (KS) probability value (Table 5). The higher 
the probability the more likely the PDFs are drawn from the same parent distribution. 
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TABLE 3 

Spearman Correlation Coefficients for Various Input Parameters and Measurement Methods 


Parameters 

R vs. a v 

a v vs. L 

R vs. L 

n = SF/F = 6a / <7 = 1 

-0.14 

0.51 

0.67 

n = SF/F = 5a/a = 2 

0.04 

0.45 

0.81 

n = 6F/F = 5a / a = 3 

0.09 

0.37 

0.87 

(edge,threshold) = (2a, 3a) 

-0.14 

0.51 

0.67 

(edge,threshold) = (2a, 4a) 

-0.09 

0.53 

0.69 

(edge,threshold) = (3cr, 4cr) 

0.06 

0.40 

0.68 

no extrapolation 

-0.12 

0.48 

0.68 

no deconvolution 

0.21 

0.51 

0.83 

no extrapolation Sz no deconvolution 

0.28 

0.48 

0.94 


TABLE 4 

Gaussian Fit Coefficients of the Distributions in Figure 16 


Parameters 

R vs 

• CTy 

Gaussian center 
a v vs. L 

R vs. L 

R vs 

Gaussian dispersion 
. a v a v vs. L R vs. L 

n = 5F/F = 5a/a = 1 

(34.39, 

14.14) 

(16.98, 

8.04) 

(31.05, 

13.04) 

(24.11, 

8.47) 

( 6.86, 

8.67) 

(16.72, 

9.27) 

n = 5F/F = 5a/a = 2 

(11.97, 

14.88) 

(16.12, 

1.40) 

(14.85, 

3.86) 

(11.62, 

10.62) 

( 8.73, 

5.97) 

( 7.02, 

5.91) 

n = 5F/F = 5a/a = 3 

(12.50, 

13.95) 

(13.93, 

1.25) 

(16.61, 

4.74) 

(14.46, 

12.73) 

(10.29, 

7.94) 

( 7.75, 

7.26) 

(edge,threshold) = (2(7, 3c) 

(34.39, 

14.14) 

(16.98, 

8.04) 

(31.05, 

13.04) 

(24.11, 

8.47) 

( 6.86, 

8.67) 

(16.72, 

9.27) 

(edge,threshold) = (2cr, 4cr) 

(33.49, 

14.42) 

(16.83, 

7.74) 

(30.00, 

13.09) 

(24.95, 

8.94) 

( 7.17, 

9.11) 

(15.70, 

9.51) 

(edge,threshold) = (3cr, 4a) 

(51.39, 

27.88) 

(38.92, 

11.43) 

(52.06, 

31.64) 

(32.53, 

14.59) 

( 9.08, 

21.79) 

(25.97, 

16.04) 

no extrapolation 

(28.00, 

34.30) 

(17.00, 

5.54) 

(25.50, 

12.18) 

(35.47, 

17.38) 

( 6.86, 

7.71) 

(22.53, 

12.72) 

no deconvolution 

(32.10, 

15.19) 

(14.96, 

7.99) 

(30.12, 

7.65) 

(14.26, 

6.91) 

( 6.54, 

8.68) 

( 9.50, 

6.04) 

no extrapolation &; no deconvolution 

(24.78, 

15.22) 

(14.99, 

5.50) 

(23.20, 

5.70) 

(14.90, 

6.88) 

( 6.55, 

7.71) 

( 8.77, 

4.66) 


TABLE 5 

Two-Dimensional Kolmogorov-Smirnov Test 


H VS. (7 V 


Parameters 

1 

2 

3 

Parameters 

(2a, 3a) 

(2a, 4a) 

( 3 cr , 4a) 

Parameters 

noex 

node 

noexdc 

i 

1.00 

- 

- 

(2a, 3a) 

1.00 

- 

- 

noex 

1.00 

- 

- 

2 

0.11 

1.00 

- 

(2a, 4a) 

1.00 

1.00 

- 

node 

0.00 

1.00 

- 

3 

0.01 

0.72 

1.00 

(3a, 4a) 

0.05 

0.04 

1.00 

noexdc 

0.08 

0.00 

1.00 






0~v 

vs. L 






Parameters 

1 

2 

3 

Parameters 

(2a, 3a) 

(2a, 4a) 

( 3 ( 7 , 4 ( 7 ) 

Parameters 

noex 

node 

noexdc 

i 

1.00 

- 

- 

(2a, 3a) 

1.00 

- 

- 

noex 

1.00 

- 

- 

2 

0.13 

1.00 

- 

(2a, 4a) 

1.00 

1.00 

- 

node 

0.00 

1.00 

- 

3 

0.15 

0.74 

1.00 

(3a, 4a) 

0.005 

0.004 

1.00 

noexdc 

0.00 

0.00 

1.00 






R 

vs. L 






Parameters 

1 

2 

3 

Parameters 

(2a, 3a) 

(2a, 4a) 

( 3 cr , 4a) 

Parameters 

noex 

node 

noexdc 

i 

1.00 

- 

- 

(2a, 3a) 

1.00 

- 

- 

noex 

1.00 

- 

- 

2 

0.07 

1.00 

- 

(2a, 4a) 

1.00 

1.00 

- 

node 

0.00 

1.00 

- 

3 

0.01 

0.74 

1.00 

(3<t, 4a) 

0.02 

0.02 

1.00 

noexdc 

0.00 

0.00 

1.00 


For the first column in Figure 16, we see that higher decomposition parameters (e.g. n = SF/F = 5a/a = 3, 
green contours) yield larger scatters than lower decomposition parameters (red and blue contours). The number of 
clouds for higher decomposition parameters is less than for lower decomposition parameters, so this larger scatter is 
not due to a larger number of data points. The larger scatter is probably due to the tendency of the program to 
combine small, neighboring clouds into bigger, merged clouds in the molecular ring island, while leaving the outer 
small islands as small clouds (Figure 1). This tendency can be seen in Figure 16, as the PDFs of higher decomposition 
parameters extend to larger radii and higher luminosities. The Gaussian centers of higher decomposition parameter 
also have larger radii and higher luminosities than those of lower decomposition parameters. In any case, none of those 
distributions shows a size-linewidth relation (—0.14 < r sp < 0.09), and they have different normalization factors than 
Larson’s relations. 

For the second column in Figure 16, the parameters (edge, threshold) = (2cr rms , 3cr rms ) almost have the same PDF 
as (edge, threshold) = (2cr rms , 4cr rms ). The KS-test yields a probability value of ~ 1. This is because we only lose a 
few clouds with peak S/N < 4cr rms . However, for parameters (edge, threshold) = (3cr rms , 4cr rms ), we lose many clouds 
because the islands only extend to connected pixels with Tb > 3cr rms (i.e. the islands get smaller), and hence the 
spread of the PDF decreases due to a smaller number of data points. None of the distributions shows a size-linewidth 
relation, and they have different normalization factors than Larson’s relations. 

For the third column in Figure 16, there is no distribution that yields the size-linewidth relation found in the Milky 
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Way. In the linewidth vs luminosity plot, the non-deconvolved distributions (red and green contours) yield larger 
velocity dispersions, as expected from equation (2) by excluding the 5u 2 /27r term, and the extrapolated distribution 
(red contour) yields higher luminosities. All distributions in the linewidth vs. luminosity plot lie below Larson’s 
relation. In the size vs luminosity plot, the non-deconvolved distribution (red contour) yields a larger size, as expected 
from equation (1) by excluding the cr 2 eam terms. All distributions in the size vs. luminosity plot lie above Larson’s 
relation. 

Overall, we conclude that there is no significant effect on our general results due to the choice of input parameters 
and measurement methods. In particular, the absence of a size-linewidth relation and the different normalization 
factors of Larson’s relations are still remain. 
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